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We present a systematic method for dealing with time dependent quantum dynamics, based on the 
quantum brachistochrone and matrix mechanics. We derive the exphcit time dependence of the 
Hamiltonian operator for a number of constrained finite systems from this formaUsm. Once this has 
been achieved we go on to calculate the wavevector as a function of time, in order to demonstrate 
the use of matrix methods witli respect to several concrete examples. Interesting results are derived 
for elliptic curves and periodic orbits on higher dimensional non-commutative geometries. 
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I. INTRODUCTION 



Recent developments in quantum mechanical theory 
have focused on the time-dependent dynamical problem, 
which proposes the use of Floquet theory in order to 
derive the time evolution operators for the dependent 
system. In our paper, we consider a new set of time- 
dependent quantum mechanical operators which define 
the physics of finite systems. Standard quantum mechan- 
ics postulates the existence of a Hamiltonian operator H 
and an initial state |^'(0)) and seeks to extrapolate to 
find our final state \'^{t)), via the von Neumann equa- 
tion. This is an initial value problem; we wish to find 
and calculate a boundary value problem that is consis- 
tent with quantum physics. This has already been cal- 
culated in the papers of Carlini et. al [l|. The principal 
methodology is to construct a time-optimal Hamiltonian 
matrix which moves us from state-to-state on the pro- 
jective manifold, in least time. These paths are geodesic 
flows; they define the surface of the state space in the 
same way that the lines of longitude on the sphere de- 
fine the globe itself. In addressing this question, we are 
naturally led to the analysis of explicitly time-dependent 
quantum systems. These dynamical systems are gener- 
ally periodic in nature for the examples we have been 
able to calculate; this periodicity allows us to use Flo- 
quet theory in order to find solutions for the problems 
in question. In considering these types of general dy- 
namical problems we find a number of time dependent 
states with constant energies; wave-vector analysis al- 
lows us to consider the geometry of the state space for 
higher-dimensional realms where it is no longer possible 
to explicitly graph the paths the state takes over time. 
Experiments have been conducted into these areas with 
the particular aim of creating these exotic qutrit states, 
as in , Q and as such it is important that we be able to 
predict best methods of state control within this regime. 
Relevant references to the geometry of qutrits and SU(3) 
operators may be found in Q.Our method differs from 



standard quantum mechanics in that we do not assume 
the form of the Hamiltonian operator; we derive it using 
the quantum brachistochrone hypothesis, being that the 
physical path is the one of least time. 



II. QUANTUM FERMAT PRINCIPLE 

This section reproduces the arguments of [l| in de- 
riving the quantum brachistochrone equation from the 
Quantum Fermat Principle, being the hypothesis that 
the state vector travels along a path of least time on the 
complex projective manifold. These results are not new 
but are included in order to outline the fundamentals 
of this methodology. Firstly, let us consider the energy 
variance on the state space. We may write this as: 



(1) 

We may define the metric of the projective state space: 

ds = AE.dt (2) 

The quantum state may be expanded in an orthonormal 
basis over a set of probabilities: 



(3) 



which then define the state space metric: 
1 



^ds^ = 1 - li-iP + d^\^)\' 



dil) 



{l~P)di') (4) 
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where P = \tp){ilj\, the standard pure state projection 
operator. For the time-dependent quantum control prob- 
lem, the action principle must be recast in such a way 
that we can apply it to find optimal transfers between 
states. Our solution, where it exists, will be a time de- 
pendent Hamiltonian that takes us time optimally from 
one state to another. We do not rule out the existence 
of systems where there are a number of extremal trajec- 
tories; however, our principal aim to search for global 
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extrema where the time is a minimum regardless of our 
particular choice of initial and final states. Writing out 
the expression for the time optimal quantum control ac- 
tion we have: 



S 



J to 



T + Cs + Cc)dt 



(5) 



ds 



AE 



Cs^i{U\i')~(i'\c^))-{U H V')+ V' H cj,)) 



(7) 



These equations are gauge invariant; we therefore make a 
particular choice of the gauge, being the reference frame 
in which the Hamiltonian is traceless: 



H ^ H 



N 



Tr{H) = H 



Tr{H) = 



The first term is unity: 



AE 
'AE 



(8) 



(9) 



(10) 



and therefore contributes a term J dt to the action, jus- 
tifying our claim of time optimality. This is the quantum 
analogue of Fermat's principle of least time. Writing the 
Euler- Lagrange equations that result from the variations 
of the action integral, we find the system of equations: 



(11) 



Defining our constraint operator: 



dH 



{{H,P}-2(h)p) 



2AE'^ 



(15) 

Taking traces and expectation values of the numerator of 
the first term in the constraint operator: 



(16) 



Tr{{H,P} -2(H) P) =0 



<{{H,P}-2(^HjP)>=0 (17) 
We derive the fundamental result: 

<i>=rr(i) = (^'1 0) + (01 ^> (18) 

Computing the anticommutator of the A-operator and 
the projector we obtain: 



P 



(19) 



We now define an associated operator to specify the 
boundary condition equations: 



G^A-(A)P 



{G{t),P{t)] =Git) 



Tr{G) = 2 < G >= 



(20) 



(21) 



(22) 



We now perform some simple algebraic manipulations 
amongst our original Euler-Lagrange equations. Left- 
multiplying (11) with {(f)\ and the conjugate of (12) with 
IV') we find: 



^) = (0|i/|V) 



(23) 



2dt^ AE^ )'^^=' 







{{H,P)-2{h)p) 



(12) 



2A£;2 



|0)(^| + |^)(0H^A,|| 

(13) 

m)^oyj (14) 



ii^i^)4(^i|( 



H 



H)i) 



AE^ 



IV) (24) 



Evaluating the time derivative of the expectation value 
of the Hamiltonian: 



dt 



<H > -- 



dH 

< —r- > 

dt 



d ~ dH 
— < H >=< — > 
dt dt 



(25) 



(26) 



This relationship allows us to differentiate under the ex- 
pectation sign. We now calculate the time derivative of 
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the overlap: 

Therefore we obtain a constant of the motion given by: 

((/) I -0) = const. = - {ipl 4)) + < A> (28) 

A similar exercise from the variational equation for the 
Hamiltonian yields: 

l'/') = ((^U) (29) 

Directly differentiating this expression with respect to 
time and substituting the constant of the motion into 
the equation we derive the fundamental law of motion: 

{i— + [G,H])\i')^0 (30) 

Some trivial algebra using Herniitian conjugates and 
right multiplication proves that the G-operator and pro- 
jection operator follow the Heisenberg equation of mo- 
tion: 

i^ = [H,G] (31) 



value, and is given by: 

foiH)^Tri^)-k^O (36) 

The other constraints which restrict the degrees of free- 
dom for the Hamiltonian of the system will often be linear 
in H: 

f,{H) = Tr{Hgj) (37) 

where the g^'s are traceless Hermitian generators of the 
unitary space. In this particular instance the Hamilto- 
nian and A-operator will obey the matrix equations: 

A = \qH + F (38) 
Tr{HF) = (39) 

F = Y.^,9, (40) 

where the Hamiltonian and associated linear constraint 
are tracefree Hermitian matrices. We take these results 
as fundamental for the time optimal evolution of quan- 
tum states and may write the equations of motion for our 
quantum dynamical system as: 



(32) 



We then use these expressions to evaluate the time 
derivative of the A-operator by direct differentiation: 



dA r - ji - d ~ 
z— = [iJ,A]+zP-<A> 



dt 



dt 



— <A>=0 
dt 



(33) 



(34) 



(35) 



It is thus proved that the G-operator, A-operator and 
projection operator all obey the Heisenberg equation of 
motion. We shall refer to this important relationship as 
the Quantum Brachistochrone Equation. This expression 
forms the foundation of applied time optimal quantum 
state control and we shall use it repeatedly throughout 
this paper. As part of our quantum control methodology 
we will consider an energetic constraint as a Lagrange 
multiplier in the action functional. This restriction holds 
the total energy used in the process to be some finite 



^-\^it))^H[t]\^it)) 



(41) 



z-{H[t] + F[t]) = H[t]F[t] - F[t]H[t] (42) 



Tr{H[t]F[t]) = 



Tr{H^[t]/2) = constant 



{G{t),P{t)}^G{t) 



(43) 



(44) 



(45) 



where G ^ {H + F) - (^H + Fj , and the derivative 

indicates an explicit differentiation of the matrices with 
respect to the time parameter. The constraint, boundary 
condition and Hamiltonian are all Hermitian matrices, 
so they will evolve unitarily via the Heisenberg equation. 
When we use co-ordinate transformations, we must use 
partial derivatives, but the distinction is made clear. Our 
formalism is explicitly time dependent, in that the ma- 
trices derived will be (in general) functions of the time 
parameter. These are more difficult dynamic problems to 
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solve than the time invariant the matrices within 

the differential equations to be solved are functions of 
time instead of constants. 



III. ISOMETRY GROUPS OF DIFFERENTIAL 
OPERATORS AND MATRICES 



IV. TIME OPTIMAL CONTROL FOR SU(2) 



For our first example, we consider the simple SU(2) 
system. In this the state space is equivalent to 

a double cover of the Bloch sphere, we expect to recover 
great circles as our solutions. Our fundamental relations 
are: 



In this paper, often we will use isometric transforma- 
tions to move between various different equivalent phys- 
ical pictures. Mainly we confine ourselves to the inter- 
action picture or the co-rotating reference frame of Flo- 
quet; as such it is quite important that we have a concrete 
method in place for understanding the meaning and ap- 
plication of isometry and similarity transforms. Firstly, 
consider the coupled matrix differential equations: 



ij^\x{t)) = {U~'HU)\x{t)} 



.dU 

'' dt 



We may then write the solution as: 

\x{t))^Q{t,0)\x{0)) 



(46) 



(47) 



(48) 



Q{t,0) = expi~i U-\s,0)His)U{s,0)ds) (49) 
Jo 



Tr{^)-k = 



Tr{Ha,) = 



Tr{H) = 



(55) 



(56) 



(57) 



We may write the expression for the Hamiltonian as the 
ansatz: 



H 



a e 

e* —a 



(58) 



Evaluating the trace of the Hamiltonian multiplied with 
the sigma-z matrix we find directly that a = 0. Our time 
dependent Hamiltonian operator is then: 



H = 



e 
e* 



(59) 



Computing the quantum brachistochrone and calculating 
the right hand side of the expression via ordinary matrix 
multiplication: 



U{s, 0) = exp{-i f V Xk{s')Pkds') (50) 



Another relevant picture we may use is the Floquet co- 
rotating reference frame. In this case, we have dynamical 
equations of motion given by: 

'Jt = + e-^^Hit)e+^') Mt)) (51) 

which has a solution that can be written in closed form 
as: 



|^(t))=A-(t, 0)1(^(0)) 



(52) 



K{t,0) = exp{-i W-\s,0)H{s)W{s,0)ds) (53) 
Jo 



W{s, 0) = exp{-~i 



dS_ 



dr) 



(54) 



i — {H + F) = HF-FH 



(60) 



we find the explicitly time dependent matrix which de- 
fines the optimal control fields to be given by: 



dt 



n e 



-n 



= 20 



-e 
e* 



(61) 



By inspection, il = const. We may write the differential 
equations for the complex control fields in the form: 



dt 



e{t) 
s*{t) 



= 2n 



-1 
1 



e{t) 
s*{t) 



Hopt (t) 



£(0)e2*^^* 
e*(0)e-2'f^* 



(62) 



(63) 



At this point the rotating wave approximation would usu- 
ally be invoked, but as the constraint is a constant of the 
motion, we may write: 

{Hopt{t)+na,)U{t,0) = U{t,0)iHopt{0)+na,) (64) 
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Using the Schrodinger equation, we find the dynamical 
equation for the unitary: 

^U{Hopt{0) + na,)-na,U (65) 

This expression Iras an explicit solution: 

U{t, 0) = exp{i^}a:,t)exp{~i{Hopt{0) + na^)t) (66) 

Some simple algebra using the identity exp{in.a) ~ 
cos 01 + in. a sin Ogives the expansion for our unitary over 
the generators of SU(2) as: 

_^ (£(2)j_^_W) sin(f7'^)[^^ sm{nt) + ay cos{nt)] 



which is of the form of the Heisenberg energy-time uncer- 
tainty principle. We may write the Hamiltonian matrix 
as: 

Hopt{t) = vi^{s\x\2Q,tax + cos2f2t(Tj,) (75) 

where = k. Hence the best we can do in terms of 
energy and time is: 

7r 

VoTmm = (76) 

In summary, we have managed to derive the energy-time 
relation for a revolving quantum state on SU(2); the arc 
that is described by our unitary time evolution operator 
has the appropriate symmetry to be considered a sinusoid 
convolved with a great circle. 



ia^ [cos{^'t) sin(m) - — cos(m) sin(17'i)] (67) 



where il' = ^/k + il?. We require a unitary that satisfies 
boundary conditions of the form: 



U{T,{))a^U\T,Q)^-a^ 



(68) 



which gives us the quantisation condition OT = m7r/2, 

TO S N. Choosing an initial state |^'(0)) = l]"^and 

v2 

a final state |1'(T)) = ^[1,-1]^, we find from the 

v2 

boundary conditions (5) the equivalent conditions on the 
control fields: 



e(0)--£*(0) 

e{T) = -e*{T) 
We may rewrite this as: 

sin(217T) = 
which gives a second quantisation condition: 
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(69) 
(70) 

(71) 
(72) 



Combining this with our first quantisation condition, we 
find the period of the wave- vector to be: 



V. UNITARY MATRICES FOR SU(2) 

For completeness, we list a number of useful unitary 
matrices that appear in the quaternionic geodesic calcu- 
lation. These are related in various ways through phase 
angles and transformations; we shall not go into this here 
as the principal part of the calculation is new work on 
SU(3). These matrices are: 



V2 



1 -1 
1 1 



Ua = 



— e 



smx —icosx 
i cos X — sin x 



1 

1 



C/5 



1 

71 



Uq = exp{—i'duz) — e 



1 

e-2»^ 



(77) 



(78) 



(79) 



(80) 



(81) 



(82) 



4fc 



(n2 - to2) 



(73) 



The left hand side of this expression is a positive number, 
and therefore n > m. The minimum time of operation is 
then: 



2Vk 



(74) 



U,Uj = U]U, = i 



(83) 



One may transform various Hamiltonians and density 
matrices to find interesting relationships between differ- 
ent arrangements of constraints: 



Fp(0 = Up{<;)H{^)UU^ 



(84) 
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" 1 e-*" 




e-^" ■ 




■ 1 1 ■ 




'1 


1 -e-"^ 




e+*" 








-1 



(85) 



VI. BILINEAR QUANTUM ALGEBRA ON 
SU(2) 



Writing the quantum brachistochrone in vector form for 
a single qubit we find: 



and hence we finally obtain the control matrices for 
the vector control system on SU(2) as: 



^2 = \/2 




1 

0-: 



0-10 



1 



(93) 



(94) 



*^(IMi)> + \im = E ('^Wl 4 l/W> (86) 

where the |j)'s are the standard basis set, and the 
matrices Aj have the property that {h{t)\Aj\f{t)) + 
{f{t)\ Aj \h{t)) = 0. Let us calculate the general SU(2) 
case as an application of the formula. In this situation we 
have system Hamiltonian and matrix constraints given 
by: 

H{t)=J2 h^ {t)c7^ and m=Y. (^)'^J- (8^) 



where the Pauli matrices have matrix representation 



1 

-1 



1 

1 



-i 

1 



Our constraints state that the Hamiltonian and F{t) are 
orthogonal, in that their dot product is zero: 



Tr{Hit)F{t)) = {h{t)\ fit)) = {f{t)\ h{t)) ^ (89) 



A, 



V2 



1 
-10 




(95) 



By explicit calculation one can show the identity 
{h{t)\ Aj \f{t)) + {f{t)\ Aj \h{t)) = holds true for aU j; 
adjoint states are defined on the space 
[ M"~^, C'^*, C'^ ] , being the complex conjugate and 
transposed row vector of the column we started with. 
Let us now apply this formulation to a well-known 
example, being the time-dependent Hamiltonian used 
by Carlini et. al in [1]. The purpose of writing the 
quantum brachistochrone in vector form is that it 
enables us to consider the physics directly, which helps 
with understanding. In this case we have Hamiltonian 
and constraint matrices as given by: 



eh{t) 
el 



Mt) 

-m 



(96) 



The vectors defined by these matrices may be written 
as: 



dt 










V2Mt) 


\hit)) = 


Shit) 


;I/W) = 















(97) 



(90) 



This may be written in the equivalent vector form as: 



V2iehe} - Sfsl) 

(91) 

Our Hamiltonian and constraint vectors obey the rela- 
tion 





' V2h, ' 










Sh 


+ 

















{fit)\A,\h{t)) + {h{t)\A,\fit))=o 



(92) 



where we have used the factor of \/2 to maintain the 
trace product relation. These are multivectors in many 
respects, having a real part and two complex 
components. The link between the matrix form and the 
vector control equation is clear; we may use either 
technique in order to calculate the time optimal 
Hamiltonian operator. 

In this situation we may write the quantum 
brachistochrone in vector form as: 



dt 



V2m 

Shit) 
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[ 72/,, 0, ] V2 
[ V2/., 0, 0] V2 
0, 0]V2 
















1 









£h 





-1 






. _ 


- 


-1 









' 













£h 


1 










Jh 








1 






' 


-1 












£h 




















' X ' 




" r " 


\h)^ 


e 


,!/> = 


n 




e* 




n* 



(107) 



{h\ = [A,e*,e] , (/| = [r,n*,n] G [M,C^] (108) 

These types of mappings must preserve inner product 
relationships between the Hamiltonian and constraint. 










£h{t) 






4W 







(98) 



dt 



as required. This system may be rewritten in the 
compact form: 



{h\D\h) ^ const. = \^ + \£\^ (99) 

(/I E \h) = E I/) = = 2Ar + e*n + eH* (100) 

U\B,\h)^-{h\B,\f) (101) 

^J^{\h) + \I)) = Y.\3){h\BAf) (102) 



' 


0" 




1 


0' 





-1 


B2=2 











1 




-1 






B. 



D = 



1 
1/2 
1/2 



(103) 











-1 




2 


1 








(104) 















" 1 















T 







(105) 








(1 


- T) 





(106) 



VII. SO(3) MATRIX ANALYSIS 

Consider a Hamiltonian matrix of the form: 



H[t] = n{t) . S = 
We form a constraint matrix 



71, (t) e{t) 

e*{t) ~n,{t) 



(109) 



The quantum brachistochrone equation reads as: 



u Ki 
K* ~2u K2 
K2 u 



(110) 



i-{H{t) + F{t)) = H{t)F{t) - F{t)H{t) (111) 

Evaluating this for our matrix system, we find the set of 
differential equations: 



u = n^=i^e* =0:H[t]=H[0] (112) 

Writing our wave-vector in standard form we find the 
relations: 



C2{t) = C2(0) 



d 

i — 






nz £ 




'c,{t)' 


dt 






£* -n. 







(113) 



(114) 



Tr{H'^/2) = const. = nl + \e\^ = (115) 
The propagator may be written in the compact form: 



U{t, 0) = i - i [0] + (cosi?t - 1) 
R 



1 

1 



(116) 
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Choosing an initial state |'0(O)) = [100] we find the 
wave- vector: 



\m) = uit,o)\m) 



cosRt — i^sinRt 


— 4e*sini?t 



(117) 



This is 27r— periodic and satisfies — |V'(^)) ~ IV'(O)) = 
I ■0(2^)^, where m is an integer. 



H{t) = R 



cos(m) 
cos(f2t) —ism{nt) 
ism{nt) 



(124) 



We now apply this Hamiltonian to a wave-vector, given 
as a superposition of energy eigenstates which may be (in 
general) time dependent. This has the expansion: 



VIII. TIME DEPENDENT ELLIPTIC 
HAMILTONIAN ON SU(3) 

In this section we shall consider three examples of SU(3) 
Hamiltonians. The physical motivation for this is to ex- 
plicitly calculate brachistochrones for some model prob- 
lems which are more general than the SU(2) models con- 
sidered in [l|. Boscain et. al has examined various sit- 
uations on SU(3) that have similar properties, see for 
example Our first Hamiltonian matrix is given by: 



H{t) = 



a{t) 
a{t) ~i/3{t) 
iP{t) 



(118) 



\^{t)) = A,{t) \+R{t)) + A,{t) \0{t)) + As{t) \-R{t)) 

(125) 



We obtain the formula for the wave-vector in^he 
standard basis \'^{t)) = [ ci{t), C2(i), C3(t) 1 as given 



by: 



ci{t) = 



sm{n't) 



Zl 



—cos{n't) — Z3sm{Q,'t) 



+cos{Qt) 



ZsCOs{il't) + Zl 



sin(rj'i) 



(126) 



with an associated constraint matrix 



Fit) 



ojiit) -t-f{t) e{t) 

l-i{t) UJ2{t) K{t) 

e*{t) K{t) uj3{t) 



= 0. (119) 



The quantum brachistochrone equation in vector form 
then reads as: 



C2(i) = [^Z2 + iR [ziCos{n't) - Z3sm{n't)]] (127) 



cos(f2'i) . , 
Zl - Z3Sm{n't) 



ism{flt) 23Cos(ri'<) + — sin(51'i) 



(128) 

where the parameters zi, Z2, z^ are given by the formulae: 



d_ 

dt 



Wl(t) 




7a 




W2(i) 


= 2 


-(7a -1- 13k) 


; (120) 


W3(i) 




13k 





Zl 




-iR 0" 




" A2(0) " 




Zl 




-iR VL 




A-(0) 


(129) 


. ^3 . 




1 




A+(0) 





d 






—j3v — [uji — cj3)q 


di 


K{t) 




-av + {uj2 - uj3)f3 



(121) 



A± = — (AiiAa) 



(130) 



d_ 
di 



a{t) 

m 



= -in 



-i 

1 



(122) ^' = VR^+n^ 

which define our initial state via: 



(131) 



d 


■ n{t) ' 







di 


. . 




aK + 



(123) 



where the complex variable e{t) — il{t) — iv{t). This 
has a solution for the time-optimal Hamiltonian which is 
given by the matrix: 



l*(0)> 



Z3 



— (fiz2 + iRzi) 
. (rizi + iRz2) 



(132) 
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In this case we were able to find a solution despite the 
explicit time dependence of the system. This example is 
particularly useful as it indicates a general method of 
solution for the quantum brachistochrone equation 
which we may use on the semisimple subgroups. 



IX. TIME DEPENDENT GEODESIC ON SU(3) 

We now take a slightly more general Hamiltonian op- 
erator and constraint than the previous example, being 
the matrices: 



where 



G{t) = (Hit) + F{t)) - m)\ (Hit) + F{t)) m)) Pit) 

(140) 

After some elementary linear algebra, we find that the 
boundary conditions are equivalent to a series of 
relations on the control fields: 



UJi — UJ2 — 



(141) 



Hit) 



eiit) 
elit) e2it) 
em 



(133) 



e2(0)=e;(0)=0 



(142) 



Fit) 



cui 

LJ2 
K* 



-(Wl + UJ2) 



(134) 



Using the quantum brachistochrone equation we ob- 
tain the set of coupled differential equations: 



(135) 



d 


■ Kit) ■ 




■ " 


di 


_K*it)_ 








d 








di 


UJ2 








(136) 



[elit)l 




-(Wl 


-uj2)ei-Ke2 




em 




(wi - 


uj2)el + K*e2 


(137) 


£2it) 




Kel - 


(wi + 2i:j2)e2 


iem\ 






+ iioi + 2uj2)e*2 . 





dt 



from which we immediately observe that Wi,W2 and n 
are constants of the motion. We now assume boundary 
conditions on the state, such that initially at time zero 
the wave- vector is in the state [ 1, 0, ] and it is in 
the state [0, 0, 1 Y after a time t = T. We then have 
boundary conditions on the projection operator as given 
by the matrices: 



P(0) 



" 1 





0" 




" 





0' 











, PiT) = 




























1 



(138) 



From [l|, we may write the boundary conditions on the 
Hamiltonian as: 



Git) ^ [Git), Pit)], Q 



<t<T 



(139) 



et(T)=ei(T) = 



(143) 



As the diagonal elements of Fit) are constant in time, 
and from the boundary conditions are initially zero, then 
we can conclude that they are zero throughout the entire 
motion. This means that the entire constraint matrix is 
constant for this situation, and we may write the time 
evolution operator as: 



Uit, 0) = expi+itF)expi-itiHiO) + F)) (144) 

Evaluating these matrices via the eigenvalue determi- 
nant, and Laufer's formulae^ we find an explicit formula 
for the state as a function of time: 



\m) 



where A 



cos(tZ\)cos(|K| t) + ^sin(|K| i)sin(<A) 



sin(iA) 



|i|Sin(|K| i)cos(iA) — ■^sin(i A)cos(|k| t) 

(145)' 



ki(o)r 



Calculating the projection operator, we obtain boundary 
conditions on the state vector: 

Pit) = \m) {m\ 



|ci(t)| ciit)cm ci(t)c^(t) 
CmC2it) \c2it)\^ C2it)cm 
cmc3it) cmc-sit) |c3(i)|' 



(146) 



After a time t — T the matrix elements of the projection 
operator satisfy: 



|ci(r)|^ = |c2(r)|2 = 0;|c3(r)r 



1 



(147) 



10 



and hence we may write the equivalent boundary condi- 
tions as: 



|c2(r)|^ = 0=(-^sin(TA))2 (148) 
thereby obtaining the condition of time quantisation: 



tion of adiabaticity nor the rotating wave approximation 
to derive our resuhs. Our method may be correctly con- 
sidered to be a competing idea, as it is technically the 
fastest way to achieve the desired transition as opposed 
to a slow, steady process as in STIRAP. 



sin(TA) = ^ T : 



ki(o)r 



(149) 



Considering the boundary conditions on the first pa- 
rameter from the wavevector. we find immediately that 
cos(T |k|) = which implies that: 



T 



(150) 



For these to both be the minimum time, they must be 
equal. From this we find the relationship: 



(n' + 1/2) 



("' + 1/2). 



ki(0)|' (151) 



The quantity on the right hand side of the equality is 
the square of a real number, and hence is positive. This 
means that the function multiplying |k| on the left hand 
side must be positive, and hence we finally derive the 
minimum time for this particular state transfer as: 



n>n' + 1/2 

TT \/37r 



(152) 



2 k 



2|ei(0)| 



After some further algebra we obtain the time-optimal 
Hamiltonian operator as given by the matrix: 



H{t) = R 



cos(fc<) 
cos(H) e~*''sin(fci) 
e+''^siii{kt) 



(153) 



where R = |£i(0)|and k = Note that as the con- 
straint was calculated to be a constant matrix, we were 
able to evaluate the time evolution operator without hav- 
ing to calculate the full time-dependent case. All that 
was required was the time dependence of the constraint, 
the initial condition of the Hamiltonian operator, and the 
boundary conditions on the state. This is a very useful 
technique, as it allows us to move to an alternative refer- 
ence frame in which for some situations the analysis may 
be more amenable to solution. Our solution is consis- 
tent with the results of Q. Much of the work that has 
been conducted on SU(3) Hamiltonians has been related 
to STIRAP 0, dl, We do not require their assump- 



TRIANGULAR MATRICES 



We now move to the more technical part of this paper, 
and derive a number of results on matrix groups. The 
reason for this is to develop a number of representations 
for angular symmetry on this elliptic group, and also to 
see if there are any general formulae which may be de- 
rived. The intent is to go from simple examples and 
steadily increase in difhculty in order to solve the main 
problems contained within. The space of complex matri- 
ces has a rich structure, and the dynamics of the semi- 
groups constitutes the behaviour of the time dependent, 
periodic state. For this reason, we use a one-at-a-time 
parameter approach, and show how each individual semi- 
group works. 

The simplest way in which a set of square matrices can be 
divided is into the values above the diagonal, the values 
below the diagonal, and the diagonal itself. Our space 
of time dependent Hamiltonian matrix operators have 
the standard Hermitian property, and this brings us to 
the consideration of triangular matrices, both upper and 
lower before we attack the main problem. Form a series 
of matrices A: 



A{a, b, c) 



1 a b 
1 c 
1 



(154) 



We now input some other co-ordinate to form the matrix 
A': 



i' = i(a',6',c') = 



1 a' b' 
1 c' 
1 



(155) 



The structure of the group is given by the anti- 
commutator and commutator, being the fundamental 
symmetric and antisymmetric operators: 



{A. a) 



1 2(a + a') (b + b') + a'c + c'a 
1 (5 + 6') + 2cc' 

1 



i(a",6",c") 



(156) 
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A, A' 



AA! - A' A : 



1 
1 




[b-b') (157) 



As we can see, these matrices form a nilpotent 
semigroup: 



" 





1 " 




" 





1 " 




" 





o" 








1 










1 













































A^{a,b,c) = 



I 2a 2b + ca 
1 2c 
1 



" 1 


a 


b ' 




" 1 


a 


b' 





1 


c 







1 


c 








1 










1 



Useful results: 





' 1 


a 


b' 




' 1 


a b + ca 


A^ =A + X = 





1 


c 


+ 





1 c 










1 







1 



(158) 



= A{2a,2b + ca,2c) (160) 



[A,X]^0 
(161) 



We can find a matrix which under the commutator 
sends upper triangular to upper triangular: 



ui 

UJ2 

0J3 



1 a{uji - UJ2) b{uji - UJ3) 

1 C{UJ2 — W3) 

1 



(167) 



An - ClA (168) 



We observe the determinant allows us to invert the 
system of equations: 







1 


a 


b 




det{A) = 





1 


c 


(159) 










1 



(169) 



Writing the differential equations in equivalent matrix 
form: 



Alt = 


1 


a 


b 




Xi 




xi + ax2 + bx2, 





1 


c 




X2 




X2 + cxa 










1 




X3 




3^3 



dx 



xi + ax2 + bx3 




dxi/dt 


X2 + CX3 




dx2 1 dt 


X3 




dx3 1 dt 



We solve the system: 



(170) 



(171) 



" 1 


a 


b ' 




' 1 


a b + ca 




" 1 


2a 


2(6 + ca) " 







1 


c 







1 c 







1 


2c 










1 







1 










1 


= x-i + naxo 

dt'' 



(162) 

The exponential matrix allows us to find the time 
dependence of the state, and hence solve the differential 
equations: 



exp(-iit) = l-iAt+ ^{-iAty 



+ l(_iit)" + ... 
n! 

(163) 



exp(-iit) = l-iAt+ —{-ity{A + X) (164) 



+ l^i-^tfA{A + X) + ^Ht)\A + XY + ... (165) 

Using the determinant matrix we find the relevant 
polynomial eigenvalue equation: 



1 - A a b 
1 - A c 
1 - A 



det(i-Ai) = [l-Xf (166) 



X2{t) = (fc + ci)e* 



(172) 



xi {t) = {^t^ + (be + ak)t + ^ e* (173) 



(174) 
(175) 



The dimensionality of the matrix group we are after, 
being the hermitian operators can be broken down into 
a sum of triangular arrays. The size of the dimension of 
the space is given by: 



1 + 2(1) = 2^-1 



2 + 2(1 + 2) = 32 - 1 



(176) 



(177) 
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3 + 2(1 + 2 + 3) = 4^-1 



n'^ ^ n + 2^ j 



(178) 



(179) 



Taking a lower triangular matrix, and evaluating the 
commutator with an upper dimensional matrix, we find 



B = 



1 
X 1 

y z 1 



(180) 



dH ~. 



dt 



dH 

(a' +/3')^ + (aA + /?/?) 



ad + /3/3 = : 



d_H 

dt 



+a 
—a i/: 


(187) 



'dt 



(188) 



H{T) = H{Q) ^ H{T) = U\T,Q)H{Q)U{T,Q) (189) 



Let 



then we find: 



— {ax + by) —bz 
— cy ax — cz bx 
ay by + cz 



a ,y = b , z 



(181) 



[AS] 



— c6* \a\ — IcI 







6*a 





ba* 
\b\' + \cf 



(182) 



Other than the diagonal entries, this matrix is close to 
the required symmetry of the elliptic matrix operators. 
The relationship between triangular matrices and 
Hermitian matrices is quite useful for understanding 
many of the problems we deal with. Any Hermitian 
matrix may be given as a sum of an upper triangular, 
its complex conjugate transpose, and a diagonal matrix. 
We now move to consideration of the particular 
Hamiltonian matrix. The Cayley-Hamilton formula for 
this matrix reads as: 

H^ = {a^+(3^)H (183) 

and we can therefore write the propagator in the form: 



exp{-ieH) = 1 + 51 (6*)^^ + 92iO)H^ 



-X a 

a —A —i/S 
i/3 -A 



-A^ + A(a2 + /32) 



(184) 



(185) 



i{ap - al3) 



1 



1 



(186) 



H{T),U{T,Q) ^ H{0),U{T,0) ^ (190) 



A periodic system has H{nT) ~ H{0) for n G N, so we 
obtain: 



H{0),U{nT,0) ^0 



U{nT, 0) = e 



-inujT -1 



(191) 



(192) 



Initially, the time evolution operator has value unity: 

;7(0,0) = i (193) 



H{t) = R 



i'^=H{t)U{t,Q) 



cos(wt) 
cos(wt) — isin(u;t) 
isin(wt) 



This operator is periodic: 



H[-(2n+l)] =i/[0] 



with initial value: 



i?[0] = 



i? 
i? 




(194) 



(195) 



(196) 



(197) 



Time commutating identities describe the actions of the 
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subgroup: 



-i 



1 



(198) 



We may directly integrate the Hamiltonian matrix to 
obtain: 



R 



sin(a;i) 



R 



cos(wr) - 1 

cos(wT) - 1 -isin(cjr) 

isin(cjr) 







sin(cjr) 
cos(wt) sin(wT) i(cos(wT) - 1) 

[ -i{cos{ujT) - 1) 

(199) 

To calculate the full time dependence of the 
wave- vector, we must exponentiate the above equation, 
which seems daunting. However, there are useful tricks 
that may be applied consistently to bring the time 
evolution operator into a more useful and physically 
intuitive format. 



XI. INVARIANTS OF SU(3) 

The invariants of a system, in a sense, define the sys- 
tem itself. Several combinations of the components of 
the wave-vector arise commonly during calculations on 
SU(3); so by considering dynamical perspectives from 
this transported reference frame we are able to see the 
crux of the physics. 



H{t) 



a 

a 

1/3 



(200) 



Redefining our system of wave-vector co-ordinates and 
constraints: 

/l = C1C3 - C1C3; /2 = C2C3 - C2C3; /s = C1C2 + C1C2; 



Im(/3) = (203) 

This means the domains of these functions are 
constrained in the complex plane. We also find in 
general that: 



dt 



dH 
~dt 



(204) 



Substituting our transformed variables into the original 
dynamical equations we obtain the system: 

f = -./4-/3/3; f = -/5-2/3N^; (205) 
dt I dt I 



dt ~ i 



/2 - /3/6; 



^ = 2«(|c2r-|ci|^)-/3/2; 



#3 , .r.df2 „ 



(206) 



(207) 



(208) 



(209) 



Defining an observable O and our algebra of states: 

'o) = UdA (210) 



1^) 



Ci(i) 
C2(i) 

C3(i) 



[cl cl c^] =(|V)*r (211) 



we derive the useful relations: 



H,0 



(212) 



/4 = C2C3 -t- C2C3; /s = cicl + c|c3; /e = C1C2 - 0^02; 

|ci|^ + |C2|^ + |C3|^ = \]a^ + P^ = const. (201) 
We immediately find the system of equations: 

Rc(/i,2) = (202) 



dt 



(213) 



After some work, we obtain the fundamental dynamics: 
dH]' 



dt 
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+ 2{aa |ci|^ + (ad + /3/3) |c2|^ + /3/3 jcal^) (214) 



Obviously this is exhaustive and sums to unity. 



dt 




(215) 



(216) 



The average of the energy, energy variance and energy 
squared may be given in terms of these functions for 
this system. 



XII. QUANTUM FRENET-SERRE FORMULA 

The Frenet-Serre formula of classical differential geom- 
etry relates the path of a particle to various properties of 
the surface it is rolling along, being the torsion and the 
curvature. This equation has a fundamental similarity to 
our Hamiltonian matrix, and as such we consider a new 
example, which is equivalent in all respects to the Frenet- 
Serre formula. By applying our method to a problem 
with known solution, we are providing a valuable check 
to our formalism. Firstly, assume the dynamical system: 



= y(v^|(l-|^) (^I)IV') (217) 
i?) = a{c[c2 + cic^) + il3{c2cl - 4c3) (218) 

i^')=«'(|ci|' + |c2|')+/3'(|c2|V|c3|') 

(219) 
(220) 



The total amount of probability is conserved: 
d 



dt 



hence we find: 



;32 « |^^|2 ^ ^2 |^^|2 ^ - C^Cs) (221) 



dt 



dt 



dt 



al3(ac2C3 + I3c{c2) = 



(222) 



+13^ =R^ aa + 13^ = (223) 



|(|C1|' + |C2|' + |C3|')=0 



a 



h + Ph = 0;224) 



To calculate the probability of any particular event, we 
use: 



P{M) = \oAt)?^\{M)? 



d 


C\ 







-iK{t) 







Cl 


I— 
dt 






iK{t) 





-iT{t) 






. '^3 . 







iT{t) 







_ ca _ 



(226) 



Writing the explicitly parameter dependent 
Hamiltonian matrix and constraint: 








-iK{t) 









a 


X 


Hit) = 


iK{t) 





-iT{t) 




a 


UJ2 


/? 







iT{t) 







. ^* 


P 


cja _ 



K,T,a,l3,uj,ef:[Q,^) 



X,X* e / : [0,oo) -> 



(227) 

(228) 

(229) 



wi + W2 + ^3 = (230) 
Evaluating the quantum brachistochrone equation: 



' dt 



H,F 



(231) 



we obtain the equivalent set of relations: 

Im(x) = const. = 77 (232) 



dr 2i 







(233) 



Hence we may write the dynamics of the curve in 
matrix form as: 





d 


'm' 




' -T]' 




■ Kit) ' 


(225) 


di 


m_ 




?7 




. Tit) _ 



(234) 
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This admits solutions: 

T{t) = ^sin(r;i) + Bcos{r]t) 



(235) 



Computing this for our Hamihonian matrix we find: 



^Pit) i/(0) = T-R [e'"*2 - e-'"*z*] (247) 



K{t) = Csin(?7t) + Ncos{r]t) 



(236) 



There is no real difference between our method and the 
Frenet-Serre method other than where (and if) you take 
the complex unit: 



dH') 
dt 



(237) 



H{t) = iA{t) 



(238) 



The solution Hamiltonian has periodicity and initial 
values given by: 



- , , ~ ,2mr ^ 

H{0)=H{ ) 

V 



In terms of the unitary this means we can write: 






~N 







N 





-B 


(239) 





B 







" 


-C 







c 





-A 


(240) 





A 








I {(^(0) \Hit)\ m)} = \ {m I [H{o).H{t)] I m 

(248) 

If we rescale the time parameter, the Hamiltonian takes 
the simpler form: 



Hit) 



-icos(t) 
icos(i) — isin(t) 
isin(t) 



Rt (249) 



The dynamical equation of state is: 

7?|^)=i?(C+|+)-C_|-)) 
which has eigenvectors given by: 



1 

71 



— icos(i) 
1 

isin(t) 



isin(t) 


icos{t) 



(250) 



(251) 



(252) 



U{T,0)H{0)U'<{T,0) = H{T) = H{0) 
Our eigenvalue equation reads as: 



-A 


-iK 





iK 


-A 


-iT 





iT 


-A 



det{H - Ai) = iK -A -iT = 
iT 

Giving us the cubic polynomial: 

- A^ + A(A'2 + T^)^0 



(241) 



(242) 



(243) 



Now. our initial Hamiltonian was constrained to finite 
energy: 



Tr{H^/2) = const. 



Hence: 



K' 



r- 



(244) 
(245) 



The solution is a circle. This fact is well known. One 
may consider various different scenarios with respect to 
the Hamiltonian-time symmetry: 



^(0) H{t) m) = (m H{o) m 



(246) 



1 

71 



icos(<) 
1 

— isin(t) 



(253) 



It is pleasing that this methodology delivers the 
standard answer to such a well-known question. The 
total curvature is a constant; the geometric figure is 
periodic. 



XIII. UNITARY TRANSFORMATIONS OF 
SU(3) 

We wish to find all the unitary transformations of 
SU(3), in order to understand the geometry of states. 
In this case, our first step is to examine the transforma- 
tions that take the time dependent Hamiltonian matrix 
from the previous worked example to a diagonal repre- 
sentation: 



HFs{t) = D{t)W{t) 



(254) 
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L ^ 



1 
0-10 




(255) 



We find: 



D{t) 



^cos(i) — ^cos(<) isnvit) 

V 2 V 2 

, V2 .V2 

;y^sin(i) --^sin(t) icos(t) 



This has the required unitary properties: 
D{t)D^{t) = D^t)D{t) = i 
with initial value given by: 

i i 

m 



(256) 



(257) 







V2 ^/2 

i 



(258) 



We may write the time evolution of this operator in 
terms of an integral equation: 



D{t) = exp{iLt)D{Q) - i / H{s)D{s)ds (259) 

Jo 

Now. consider the eigenvectors of our elliptic 
Hamiltonian matrix: 



cosiji 
cos(f> 
isinri 





-isintZ 




1 


COS(p 


1 


COS(/) 


1 


1 


72 


ism4> 


'72 








COS(f) 







— sin0 




" 0' 


COS0 





— isinc/) 






















icos(j) 








(260) 



(261) 



cos(/; 
COS0 
isinri 








1 


— COS(/) 


-1 


— COS0 




1 




1 


71 


— isin0 


"71 


— isin0 



(262) 

Writing a matrix with these eigenvectors as columns, we 
obtain: 



J(0) 



1 

71 



1 



C0S(7 



1 
71 



V2 



COSfi 



a/2 



-smi^ 




71^ 



V2 



smffl icosG 



(263) 



This matrix is unitary: 

j(</))jt(</>) = jt(0)j(0) = i 



J{(f)LJ^{4>) 



cos^ 
cost/) — isin^ 
isin(/) 



(264) 



(265) 



r 1 



J(0) = 



1 







^/2 ^/2 

V2 y/2 
i 



(266) 



and has initial value given by a Hadamard gate with a 
■J rotation on the ancilla. Finally, using the completely 
complex elliptic Hamiltonian 



(267) 



Eigenvectors: 



& 



1 

71 



V2 





' ei 


" 




H = 


el 


£2 






e* 







C0Siy9 




ie~ 


*^siniy9 


1 


& 











COSip 



(268) 



—cosip 

1 



, Q = const. G [0, 2tt) 



7i;"^ 


--Leos^ 


ie '^sin(^ 


^ 71 


71 





-^e-sin^ 


* in ■ 

=e '^smcp 

V2 


COS(^ 



(269) 



0(0) 



7= 

a/2 a/2 
1 I 

V2 ^2 
1 



(270) 



The initial value of this operator is a Hadamard matrix 
with an ancilla bit set to 1. It is unitary: 



Qt(^)g(^) = g(^)gt(^)^l 



g(^)iQt(^) = 





C0Siy9 



COS(p 





-sm(y9 



ie^^svatf 



(271) 



(272) 
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We may summarise these formulae on unitaries into the 
compact identities: 

S{k[t];t)^ [ k(s)Lk'<{s)ds (273) 
Jo 







— 2isin6'cos6' 
-isin6'cos6'(cos6' — 1) 



-2isin6'cos6' i{cos6 — l)sin6'cos6' 

,2/1 



sm 




(281) 



§-0^§Lk^^kL'-^. (274) 
6A dA SA 



H, k Lk^ = kL H, k^ (275) 



kLk^^ = k l| k^ (276) 

Our propagator for the fuh Hamiltonian may be written 
in the compact form: 



Uie,0) exp{-i / dsH{s)) 
Jq 



1 + (c — l)c^ ~isc — i(c— l)sc 
—isc c —s^ 

i(c-l)sc l + (c-l)s2 



c = cos6', s = sin6' 



(277) 



L,J{4>) 



kQiifi) 



V2 





V2 
sinri 



-\/2cos(/ 




-smq 





(282) 





-V2 



^2 



e'^siniyS 



\/2coaip ie *^siniy9 





-^e*^siniy9 



(283) 



i-\/2cosx isinx 
-V2 

-^sinx -TT^sinx 



V2' 



(284) 



Hence our fundamental rotations form compact 
semi- groups with the operator L. One may define a 
sequence of polynomials via the time evolution 
operator, viz.: 



We have eigenreflections: 



q{cos6) 



(285) 



Lf''{0,O)U{e,O) = C/(6i,0)C/t(6i,0) = i 



(278) 



\vk) {Vk\ 



(286) 



U{e, 0) = i - ism9H + {cos9 - l)H^ 



(279) 



One final useful unitary transformation is the equivalent 
NOT-gate, which can be used to flip two state labels: 



Ah 



1 



-cos^ 
isin^cos^ 







-cos^ isin^cos^ 

1 isin^ 
-isin^ 1 + cos^^ 



— isin^cos^ 

1 



isin^cos^ 1 + cos^^ 



(287) 



(288) 



N 



' 


ei 


" 




" 




" 




' 





1 " 











£2 





£i 


■,N = 





1 








p* 

^2 










£1 







1 









(280) 



Writing |fj)for the j-th column of U, it is simple to 
show that {vj \vk) = Sjk- Any rotation in this space is 
composed of these fundamental operators . We now 
compute commutators of our fundamental 
transformations with the matrix L: 



1 + sin^S^ cos^ isin^cos^ 

cos^ 1 — isin^ 

— isin^cos^ isin^ 1 + cos^^ 



(289) 



These eigenflections obey the periodicity condition: 

Af2(-C) -Mi(0+M3(0 (290) 
which is directly related to the equation of state: 
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|3)(t+j) - = V2 |2>(4) (291) 6q(s) = -2(5^2 _^ ;l20s1° + 5016s** + 86527s^ 



XIV. ELLIPTIC POLYNOMIALS 



550413s* + 895923s2 + 396900) (303) 



Using the formula for the propagator, and making the 
substitution z = cos0 we obtain the associated 
z-polynomials: 



(292) 



Pi.™(z) =q(z) 

We display several of these functions: 

q{z) = 4z^ - - 8^^ + + 42^ -62^ + 1 (293) 



s(s2 + l2)(s2 + 22) X • • • X (s2 + 7^) = s n (s2 + fc2) 



fc=i 



(304) 



6p(s) ^ -(s** + 29s'5 + 2083-* + 306s2 - 144) (305) 



Pi,i(z) = 22-*- z-^ -3z2 + 1 =39(z) (294) 



P3 3 (z) = _2z4 + + z^ - z = P(z) (295) 



rp(s) = rp(s) = sj|(s2 + fc2) (306) 



fc=i 



For a small circle centred on the origin: 



rp{z) 4 



(307) 



qiz) + PsA^) ^ Qiz) 



(296) 



g(z) = 4z^ - 4z^ - 8z^ + 6z* + 5z^ -5z^ - z + 1 (297) 
We may Laplace z-transform: 

L[F{z)]^ Lf[s]^ F[z{0)]e-"'^de (298) 
Jo 

These transformations enable us to move between the 
different representations in this space. They are 
analogous to the matrices used earlier to transform to 
the eigenbasis. These transformed functionals are 
rational polynomials with various interesting properties: 



Mf)rf,. 1 S'-^,dz^-'- (308) 
TQiz) 2 J rp{z) 4 

r r 

We may further Fourier transform our functions to find 
e.g.: 



Fw{k) = / e-''''Lw{s)ds 
Jo 



Fp{k) = / e-'^'Lp{s)ds 
Jo 



Us] = 



Lq[s 



t r„i his) 



Lp[s] 



(299) 



(300) 



H[k] ^ { I 



0;a: < 
a- > 



bg{s) = s^ - 125s^ + 24s* + 1923^ - 960s2 - 2880s + 20160 

(301) 



(309) 



(310) 



Fp(fc) = — {u[^k]H[k] - u[k]H[^k]) (311) 



u[k] = e*'^ - 6^'= + 26.2'= + e'^ + 1 (312) 



(313) 



The other transformation we require is given by the 
Mellin integral: 



r,(s) = s« 



(302) 



M[a] = / z"-V(z)dz 
Jo 



(314) 
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One immediately obtains the p-adic sequences: 

MFLp[a] = M[Fp{z)l MFLQ[a] = MiFgiz)] 

MpLpla] = M[Fp{z)], 



order must not be interchanged except in unique 
circumstances. These functions, and their related theta 
transforms are part of a general group of elliptic 
symmetry operators. One may plot these simply on a 
polar plot diagram. 



AWH = ^r(a)(ij 



— — 

3a + 2" + 



1+3 



i7rS(a) 



XV. SEMIGROUP TRANSFORMATIONS 



(315) Consider the unitary matrices ^Q,D, j|. We may 

write these matrices as the outer product of unit 
eigenvectors with column vectors: 

i7rS(a) 



M[Fq{z)] ^ -Via) 



1 2 



1 22 

5" + 2" 



(316) 
(317) 



Q = 



qi 92 93 



[\qi) 192) k3>] (324) 



Mellin integrating the original Laplace transform, we 
find for lal < 1: 



MLp[a] = -- ^ ' 



MLQ[a] 



7r(a — 1) 



Y{a) 



96 . / n(a - 1) 
sm ' 



TT X{a) 



8 . fn(a~ 1) 
sm ' 



MLp[a] = - 



X{a) = f-3 + 2" - 3"-i + 22("-i; 



(318) 



(319) 



(320) 



(321) 



J = 



Jl 32 J3 



\h) \h) 



(325) 



D 



di d2 dz 



[Ml 



M3 



(326) 



l.m 



^^„,{s). We first 



Writing |^/(s))„ = [A{s) 
consider the various transformations on the unitaries: 

L"-(a)|a,((7)) = |a,(a + a)) (327) 

Finding the matrix operators: 



coscr sincr 

1 
—sine coscr 



coscr —sincr 

1 
+sincr coscr 



Y{a) = (3 + 33 X 2" - 3" + 3 X S^""^) + 6" - 3 x 7'-°'-^^) 

(322) 



(328) 



Zia) 



1 + 2" - 3° 



22(a- 



(323) 



The other rational function has a pole at the origin. In 
this description of the space the Fourier, Laplace, 
Mellin integrals and cosine transform play a 
complementary role to the unitary matrices used earlier. 
The order of integration is compositional, in that we are 
fully integrating over the variables at each time and the 



coscr ie sincr 
1 
ie'^sincr coscr 



coscr ie sincr 
1 
-ie**sincr coscr 



(329) 
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XVII. ELLIPTIC WAVE EQUATIONS 



coscr isincr 

1 
zsiriCT cosCT 



coscr —isincr 

1 
-isina cosct 



coscr —isincr 

1 (330) 
isincr coscr 



XVI. FERMAT PRINCIPLE 

Take the action function given as the classical Fermat 
principle: 



Writing \Ai{s))^ 



A{s) 



l.m 



we now 



consider various operators on the respective 
wave-functions. It is important to note that the variable 
may properly be held to be a periodic complex variable, 
and the function extended meromorphically to the 
complex space. We will work only on the section where 
the wave is continuous, invertible and free of poles, so 
we may multiply and divide by trigonometric functions. 
In this section, as we are dealing with changes of 
variables, we must use partial differential operators to 
distinguish them from the strict time derivative treated 
elsewhere in the paper. Let us examine a particular 
co-ordinate transform, to show how we may write down 
the wave equations. We have: 



S{ya,ya\t')= n{ya)\/yldt 
Jo 

The Euler-Lagrange equations then read as: 



/Tj On _ d_ ( yan{ya) \ ^ 
" dya ds \ . /7i^ 



(331) 



(332) 



We consider the particular case of one co-ordinate: 

n{ya) = n{y) (333) 

and transform the co-ordinate as y = ia sin The dif- 
ferential equation then reads as: 



dn d4> dn{y) 
dy dt ds 



= 



(334) 



We may therefore write, after taking the appropriate 
limit, that: 



dn 



jdcj) dy 



(335) 



dn 



In general, — — ^ 0, so we therefore have the parametric 
ox 

identity: 



(336) 



dt ^a'^ + yi dy 



Applying this twice, we find the elliptic differential equa- 
tion: 



9^ 



1 



d 



1 



^ y2 dy + jy2 dy 



(337) 



[i?(x)]3,i = ^sin(x) = ri^i = r (339) 

where we drop the subscripts for convenience. Directly 
differentiating: 



dr 



d^ = 71'°^^ 



Squaring the derivative: 



{dT)'^-{-+r^)dx' 

We find after inverting the fraction and taking the square 
root: 



dx ^dx 
dr dr 



(340) 



2+" 



Partial differentiating a function twice with respect to 
our parameter: 



dxY d^^ 54- d fdx 



dr"^ \dT ) dx^ ' dx dr \dT ^ 
Substitution of the appropriate values brings: 



(341) 



2 ^ ' 



dx^ 



■ dx 



(342) 



, 2 ^.d''(t> 92 



y d(j) 



dy^ v/a2 + y2 dy 



(338) We can then find the differential-integral equations for 
the wave-function. In general it takes the form: 
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M = [0] + <3 r (343) 



If z = r cos 9, f{z) = + iv{z), we find the 
symmetrical equations: 

5;: = -cot^^ (354) 



e,^^. e {l,-l,*,-z};/35 e {1,-1}; = — (345) 

i^2''^''°) e {(cotj/, cos^i/), (tani/, sinV)} (346) 
Our entire wavefunction is then: 



(347) 



(348) 



These formulae are generated by a group of partial 
differential equations of type: 



92 f 



9* 



a2 9x 



z = fix), \z\ < a, fix + moj) = fix) 



(349) 



(350) 



and we take the angular variable so as to avoid any 
discontinuities in the relevant trigonometric ratios other 
than at the extremum, where we are guaranteed 
analytic continuity by the hyperbolic sine-sine 
relationship. They are on the perimeter of the 
wavefunction; we are guaranteed continuity and 
unitarity by virtue of our matrix construction. One can 
label the functions by their Fl'^'^ ^ ^ and . If we 
consider these functions as proper analytic functions in 
the complex plane, it is straightforward to derive 
various differential systems, analogous to the 
Cauchy-Riemann equations. For example, if z = re*^, 
f{z) = u{z) + iv{z), we obtain the cylindrical polar set: 



du 
dr 

dv 
dr 



1 dv 
rde 

1 du 
'rde 



5V 19/ 1 

Qj.2 J. Qj, j,2 Q02 



(351) 



(352) 



(353) 



dv tO^^ 
dr do 



(355) 



If z = ircosO, f{z) = u(z) + iv{z), we find polar elliptic 
equations: 



du 
dr 



dv 



dv 



^/z^~fr^ d9 



du 



dr ^2 de 



(356) 



(357) 



Differentiating the correct parameterisation, 
analytically continuing it to the complex plane and 
using separation of variables allows us to solve for the 
ground state of the system, which can be quite difficult. 
The radial part of the differential equation is generally 
the circular Bessel function: 



d'^R{r) , ldR{r) ^ a 
r 



+ II^(^) (358) 



c?r2 r dr 
The angular part of the particular solution is then: 

'^{r,e) = R{r)Y{e) (359) 



+ ^i(^)^ + ^FMYiO) = (360) 



XVIII. INTEGRAL RELATIONS 

Expanding the Hamiltonian matrix in the spectral rep- 
resentation, we find: 



(361) 



In this particular situation, the eigenenergies are equally 
separated from zero, we may therefore write: 

# = i?(|0+)(0+|-|0_)(0_|) (362) 

An alternative way of writing this Hamiltonian is: 

H = Y,^,,{t)Q,, ^ e:^{t)Ql (363) 



^,3 



(364) 
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Using the dynamical equations of motion, 
. d 



dt 



E 



1 
-1 



(365) 



Let us now expand the wave- vector in the eigenbasis 
d 



dt 



(366) 



Dropping the sums for convenience and assuming Ein- 
stein summation convention, we evaluate the differential 
equation: 



(367) 



^cAt)J^\m) = {E,^^^)\m) 
Consider a path integral on our state space: 

[Kjk{TM^ jj^ ^M-i J dsC[E,,c,\s])D[x^is)] 

(368) 

We may write the Lagrangian density as the formula: 

C = ^E,cj{s)5,k + ic,{s) (fc(s)l H{s) \j{s)) (369) 

Up to isomorphism, this fulfils an equivalent role to the 
time evolution operator. As we have managed to already 
find a number of differential relationships, the necessary 
question is to enquire as to whether there are a series of 
equivalent integral identities that we may use to evaluate 
these particular groups of elliptic differential equations. 
Consider a simple toy problem of a wave-packet interact- 
ing with a time dependent oscillating potential: 



ft2 a2*(a;,i) 
2m dx^ 



VQCQ'&ujt.^ix.t) (370) 



Using natural units and rescaling the time, we may write 
this as a neat operator equation: 

*^^^ = (^o + ^/(i))*(^,<) (371) 

Now. if we cosine-transform the time co-ordinate, our 
equations of motion transform: 



z = cosa;< 



t = — cos ^(z) 



Hi{t) ^ Hi{z) = Voz 



(372) 

(373) 

(374) 
(375) 



d dz d I d 

Writing our original wave-equation in the transformed 
frame, we find: 



zo; vl 



dz 



2 



+ Vqz^(x,z) (377) 



Using the method of separation of variables, we write an 
ansatz for our wavefunction: 

^(x,z) = (\){x)ii{z) (378) 

Finding the equivalent set of equations: 

1 d^(b{x\ 



E(\)(x) 



2 dx^ 



dz 



(379) 



(380) 



The first of these equations is readily solvable using com- 
plex exponentials, however, the second presents us with 
some difficulties. Firstly, we rescale the variables to place 
the differential equation into dimensionless form, finding: 



dz y/T^"^^ ' 

Directly differentiating: 



(1 - z^) 



(Pi 
dz^ 



{p{a, z) 



(381) 



(382) 



where the p's and q's are polynomial functions. This is 
one of our basic elliptic differential equations. Note the 
useful integral formula: 



Vb 

/seed + tan z)dz = / 



dy 



i{yb-ya) (383) 



where y = secz. To consider further transformations, 
consider first one of the polynomials from the previous 
section: 

hi (z) = + 29z^ 4- 208z'' + SOGz^ - 144 (384) 

It is possible to completely specify this polynomial by its 
root set, viz.: 



6i(z) = 0, z € {±q;i, ±i/3i, ±i5i, ±171} 



(385) 



Moving sequentially through the other relevant polyno- 
mials, we find: 

62(z) = z^-125z5 + 25z^ + 192z3-960z2-2880z-f 20160 

(386) 
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b2iz) = 0, z e {a2,/?2, -72,^1,^^,^2,22} 



(387) 



This point set has the neat property of being a pentagram 
in the complex plane with three non-symmetrical points 
on the real axis, one of which is the primary vertice of 
the pentagon. 



63(2) = z*^ + 120z'"^ + 5016z4 + 865272^ 



550413z2 + 896923Z + 396900 



(388) 



This polynomial has its root set confined to the negative 
real axis: 

63(2) = 0, z e {-as, -(33, -73, -S3, -63, -Ca} (389) 

Finally, consider the polynomial 64(2) = 6i(\/z).As the 
original polynomial is even in z, the new square-root 
transformed polynomial is of the form: 

bi{z) = z^ + 29z3 + 208^2 + 306z - 144 (390) 

which has a root set: 

64(2) = 0,ze {a4,-/34,-74,-M (391) 

By examining the cosine and exponentially transformed 
polynomials, bj{cos0) & bj{cxpi9), and real and imag- 
inary parts thereof, it is a simple exercise to develop a 
series of polar representations for these polynomial func- 
tions. The complex analytic behaviour is generated by 
these or similar sets of roulettes. They are related to a 
number of interesting integral relations. We begin with 
a simple example: 



271- 



hi{u)du . „, . r—T, 



(392) 



The other polynomials have similar identities with re- 
spect to this measure: 



bi{u)du _ 1233l7r 
1 " 128 



(393) 



37r 



6i(w)Vl - u'^du = bn+i\/9^^^~~l{c+{9n^-l)P{TT)) 

(394) 

3 

(395) 



P(^)^^C„^2"+i 



27r 



bi{u)x/l-u^du = A:7r+V47r2 - 1{1+{4:Tt^-1)Y{tt)) 

(396) 



71=0 



{Dn}C{Cn},Dn^OmodCj 



(397) 



(398) 



tt/2 



bi{u)\/l - v?du = gix+i^Tx'^ - 4(/+ (tt^ -4)r(7r)) 

(399) 

3 

ri=0 



{M„} C {C„} 



(401) 



bl{u)^Jl - u^du = - Z2(^C„7^"+^)-msin-^(0 



(402) 

r ^^jM^l^ip{a)V^^^+Xa + ism-\a) (403) 
7o V 1 — 

Off the unit circle, these integrals lose analyticity unless 
we modify the weighting function. One way to do this is 
by defining a co-distribution: 



+1 r(^ + i) 

-1 



(404) 



We then form marginal distributions: 



J^Pkis)W{s,a)ds 
In ^(^' ^)ds 



(405) 



6i,w(a) 



= -3 



p(a) 



ll{a + 2{n+l) + l) q{a) 
Using partial fractions to expand the denominator: 
1 _ 1 ^ 1 1 



(406) 



n(a + 2(?i + l) + l) 48 a + 3 a + 9 



1 



16 a + 7 a + 5 



(407) 



p{a) = 48a'* + 1050a^ -I- Ib'd^a^ + 17653a - 1214 (408) 



lim 61,1^(0!) = -144 



(409) 
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This is not the only co-distribution, we may form an 
equivalent functional: 

W'{a, s) = (1 - a2)-.V2 = W{-s, a) (410) 

These integral relationships are then neatly expressible 
in terms of hypergeometric functions: 



1 2 

daW{—s, a) = a < 2 ^ 

2 2 



(411) 



r 1/2 ] 

daW{-s,a)bi^w{a) = 2^ < f !1 + 1 f Q^" 

ri=l L 2 2 J 



(412) 



The point symmetry group of p{a) is a kite, all the zeroes 
of q{a) lie along the negative real axis. 



XIX. CHEBYSHEV POLYNOMIALS 

Consider a unit triangle in the complex plane. Using 
De-Moivre's theorem: 



(e'^)" = {cos9 + ismt 



= cos m6 + i sin mfc 



(413) 



We may write an alternative transformation, viz. 



z = X + iy = X + i\J \ — x^ (414) 
The left-hand side of De-Moivre's theorem then reads as: 



(x + i\/ 1 — x'^y^ = cosmO + isiimiO (415) 

We define the Chebyshev polynomial of the first kind: 

T„i(cos6') = cosm^ (416) 

Using the double angle formula, it is a simple exercise to 
derive the first recursion relation: 

T,nix) = xT,n-l{x) + (X^ - l)U,n-l{x) (417) 

where we have defined the Chebyshev polynomial of the 
second kind: 



Ujnicos9) 



sin m6 



sm ( 



with recursion relation: 



(418) 



(419) 



Directly differentiating the recursion relations, we obtain 
the system: 



dx 



ml/mix) 



(420) 



(l-x^) 



dUm 

dx 



il-x') 



d^T 
dx^ 



m{xUm - mTrn) (421) 



{mxUm - m^T^) (422) 



[l-x^y- 



d^T 
dx^ 



dTm 
dx 







(423) 



This differential equation is related to the elliptic wave 
equations used in the previous sections. This comes 
about due to our particular break-down of the wave- 
function in terms of roots of unity on the unit circle. 



XX. CIRCULAR CO-ORDINATES 



Consider the standard cylindrical co-ordinate system 
in two dimensions, i.e. a circular system: 



cos 6 
sin 9 



(424) 



Directly differentiating to find the velocity vector: 
dx 



dt 



rir + rOig 

The polar unit vectors have representation 
^^(^)=(srnl)'^«W = 



(425) 



— sin 6* 

008 6* 



er.B{9 + a) = L{a)erM{9) 



L{a) 



cos a — sm a 
sin a cos a 



Writing out the arc-length for the system: 

ds^ = dr"^ + r'^d9'^ 
we derive the gradient for the wave-function: 

^ dr r 89 
Deriving the Laplacian in circular co-ordinates: 

V^* = (V-V)5' 



Qj,2 J, Qj, j,2 QQ2 



(426) 



(427) 



(428) 



(429) 



(430) 



(431) 



(432) 
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The time-independent wave equation for a free particle 
on the circle is given by 



H^{r, 9) = -V^1'(r, 9) = ^^'(r, ( 



(433) 



Using the method of separation of variables, we obtain 
the equivalent system: 



■^{r,e) = R{r)Y{9) 



(PY{9) 
d9^ 



(434) 



(435) 



'^^'-'-^HK-Wr)^ (436) 
ar^ r ar 

These equations have the well-known solutions of com- 
plex exponentials and Bessel functions. The angular part 
obeys periodic boundary conditions, and can be sepa- 
rated into an odd and even part, being the sine and co- 
sine: 



Y{9) = Ae"^" + Be 



Y±{e) = ±Y±[ 



Y{2n) = y(o) 



ikS 



(437) 



(438) 



(439) 



finding a particular quantisation for the system as given 
by: 



k = ±7 



(440) 



Rewriting the radial wave equation in natural units ^ = 1, 
we find: 



r — — — \-r ' 



[n" - r^)R.^ = (441) 



dr^ dr 
Our solution is then the Bessel function: 

Rnir) = Jn{r) (442) 

This has an integral expansion and inner product: 

Jn{r) ^ r e-*("'^~''''''^'^)d(/> (443) 
27r ,1^ 



_ sinijiim - n)) _ S{m - n)S„ 



2'K^{m — n) 



27r2 



(444) 



XXI. SPIN WAVES IN CRYSTALS 

Consider a spin wave that moves in a one dimensional 
crystal, perhaps along a line. This is the simplest of the 
continuous degrees of freedom that can be analysed; we 
examine it purely to explain the context of the Bessel 
function that arose in the previous section, to illustrate 
the direct physics of the phenomena in quantum systems 
that this paper is concerned with. We take the Hamilto- 



H\n) ^ -A{\n + l) + \n-l) -2\n)) (445) 
Defining a matrix representation: 

Hra,n = {n\ H |m) (446) 



(447) 



We solve this differential equation with the ansatz: 

Cmit) = a™e-^^* (448) 



EUn = A{2an — a„_i — Un+l) 



AE{p) ^ E- Eq = 2cos(p) 



(449) 



(450) 



This eigenvalue equation has a solution that can be writ- 
ten in integral form as a Green's function: 

1 

K{Aq,t\0,0) ^ — / exp{-i{pAq -2t cos p)). dp 
27r J_„ 

(451) 

K{Aq, t\0) = {-i)^UAg{2t) (452) 
This is the Bessel function, as required. 



XXII. TIME DEPENDENT OSCILLATOR 
REVISITED 



Take the wave equation in one space dimension, cou- 
pled to a time dependent oscillator: 



ih 



2m„ dx^ 



Vo cosojt.^' 



(453) 



Using the cosine transform, and the method of separation 
of variables: 



(454) 
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^ 1 _i, , 

Z = cos Ujt ^ t = — COS (Z) 
LO 

we find the differential equation of state: 
.dp (1 — az) 



which has an exphcit solution: 



(455) 



(456) 



i(ln^(z) -ln/3o) = - cos"^ [z) + ay/ 1 - 2^ (457) 

This is mid-way between a Bessel function and a Cheby- 
shev polynomial. Writing x = "^^j ^ = cos xwe may show 
a neat set of transformations: 



9* 



-1 9^- -1 94- 



az Vl - 2^ sin X dx 



(458) 
(459) 



Differentiating these expressions directly, we arrive at two 
equivalent differential equations: 



d 



„ , - v/l-z2 — (^1-22*) (460) 
ox oz 



lfT = ^^i^7r^ (461) 
az^ sm X Ox sm x ox 

Expanding these, we find the pair of second order differ- 
ential equations, related by the transform pair: 



dx^ dz"^ dz 



(1-z^) 



cot 

ax^ ^ dx 



a* 



a^f 



oz ax 



(462) 



(463) 



(464) 



Taking the equation of state and directly differentiating 
we find a second order differential equation: 

n^d^B , ,n I z(l — az), 

i{\ - z2 = /3 -z 1 - azf - aVl-22 + ^7== 
az-^ VI - z^ 

(465) 

This may be rewritten in the equivalent form: 

(1-^ )tT + (^ i-z)^ + {l-azfp = (466) 

az^ 1 — az az 



/3(z) = /3oea:p(« cos ^(z))ea;p(— ia-\/l — z^) (467) 



XXIII. DISCRETE FOURIER TRANSFORM 
FOR SU(3) 

The discrete Fourier transform (DFT) on a qutrit has 
subtleties that are not apparent in the even dimensional 
case, as there are several transformations which have sim- 
ilar properties to the expected form. As the structure of 
the entire unitary space can be broken down into three 
sets of orthogonal rotations and two diagonal operators, 
defining exactly what is meant by 'odd' and 'even' is not 
quite so simple. This is due to the underlying group 
structure, as we may define a third set of properties, 
namely neutral. We would expect to have three sub- 
groups which map in a similar way to the DFT on even 
dimensions, and one total unitary; this is indeed the case. 
To show that this is true, we first generate a cube root 
of unity: 



(468) 



with complex conjugate: 



z* = l(-l-zV3) 

This has unit modulus and is a group under {*,2 }: 
|z|^ = l 

z'' ^z\z^{z*Y 
Defining our first DFT matrix: 
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71 
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1 " 


1 

^71 
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1 ' 


1 
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z 
z* 


z* 
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1 
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z 

Z2 


z2 

z4 



This matrix is unitary: 

RR^ = Ft^R ^ i 



(469) 

(470) 
(471) 

(472) 

(473) 



and is a fourth- root of unity = 1. One useful relation- 
ship is: 



RF R = Rlx 



1 
1 
1 



(474) 



which is the equivalent C-NOT for qutrits. We may write 
our complex unit as the cube root of unity: 



.3 = 1 



with equivalent Euler representation: 

2ni 

z = exp( — ) 



(475) 



(476) 
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Let us now consider a particular decomposition of the 
Hamiltonian: 
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^Ho + Hq 
(477) 

We define an odd- valued split DFT via: 
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1 1 f 
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with function values: 



A: = e'^fc-i = fc* (479) 
This takes our Hamiltonian to the operator: 



H 



cji 1 - fc-i 1 - k-^ 
l-k UJ2 1 - 
1-k^ 1-k^ ujz 



(480) 



Writing out the formulae for the trigonometric polyno- 
mials: 



cos 6'(4 cos^ e-i) + i sin 61(4 sin^ 61-3) 



k'^ = 2003^6* - 1 + i2sin9cos9 



il) 



(fc")* 



(482) 
(483) 



Note that this takes the Hamiltonian to the off-diagonal 
non-comutative submanifold; consequently we expect to 
find two equivalent different split DFTs which take the 
Hamiltonian to the same type of submanifolds. This type 
of transformation will not be unitary or invertible, it is 
more along the same lines as making a projection onto 
a co-ordinate axis, keeping a subset of the co-ordinates 
and discarding the rest of the set. Writing out the trans- 
formation law on the Hamiltonian matrix: 



Hw = W 



1 
0-10 




(484) 



We now substitute the even-valued split DFT into this 
relation: 



W 



1 1 1 

1 w w'^ 



(485) 



and find the explicitly transformed matrix: 

-2 ■ 



Hw — 



1 

I — w 



- w 




-1 1 
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1 



- w 




(486) 

w = e*'^ (487) 

This matrix has two independent component func- 
tions, thereby satisfying the physical conditions of the 
SU(3)/SU(2) submanifold. We may write this in the 
form: 

Hw ^ A + A^ Xcosx + y sinx + K + A"^ (488) 
Consider the operators X,Y: 
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(489) 



They form a non-commutative subgroup that is similar 
to SU(2): 



-2* 



1 


0-1 



2iZ 



[X,Z] = -iY 
[Y, Z] = iX 



(490) 

(491) 
(492) 



Up to a constant matrix, this Hamiltonian is equivalent 
to the S0(3) subalgebra of SU(3). The constant matrix 
is the set of upper and lower projectors, with representa- 
tion: 



A = Ai + A2 



A. 



(493) 



ii = |l)(2|,i2 = |l)(3|,i3-|2)(3| (494) 



(495) 



This matrix is not unique, one may define another ma- 
trix, given by: 



J = 



111 

1 j f 
1 .f f 



(496) 
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Transforming the Hamiltonian we find: 



Hj = J 







1 
0-10 




Jt 



Hi 



1 - i 1 - j-^ 

l-f I- J 



(497) 
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(507) 



This is the smallest non-commutative group. Note that 
these are a group of rotoflections, as opposed to proper 
reflections or rotation matrices. Using our Hamiltonian 
operators, we may deflne a closely related series of oper- 
ators using our initial value equations: 



which means our mapped Hamiltonian is an equivalent 
two-parameter subgroup. Writing this as a congruence 
relation: 



Hq ¥ Hj ^ H 



w 



(499) 



The functions involved have decompositions that may be 
readily calculated with De-Moivre's theorem: 

1 - j2 = 1 - e^''^ = 1 - cos 261 - i sin 29 (500) 



1 — cos 29 — i sin 26' = 2(cos2 6l-l)-i2sin6'cos6' (501) 



1 - r 



2(sin^ 



i2 sin 6 cos 6 



(502) 



These functions have the necessary symmetry under con- 
jugation to give us a Hermitian matrix: 



irr 



(503) 



Q(0) = 



1 







V2 V2 

V2 ^2 
1 



(508) 



This is one of our equivalent Hadamard gates, which con- 
stitutes a pure rotation in 3-D space. 



0^(0) 



0-10 
1 
1 



(509) 



We then have a reflection in a plane of symmetry, which 
is a signed permutation matrix. 



^(0) 



>/2 \/2 ^ 

V2 V2 ^ 
i 



(510) 



i-r' = (i-j')* 



(504) 



None of the split-transforms J, Qor W are unitary or or- 
thogonal; the only matrix which has the required prop- 
erty to be a DFT in its own right is R. 



Because we are working on higher dimensions, there is a 
greater number of available operators, which we need in 
order to close the algebra and conduct meaningful calcu- 
lations. With qubits all that is required is a phase gate, 
a NOT gate, a Hadamard matrix and a XOR gate. 



i-1 i+1 
+ 1 i + 1 
0-2 



(511) 



XXIV. DIHEDRAL GROUP ON THE UNIT 
TRIANGLE 



The permutations one may apply to a triangle with 
uniquely labelled vertices generate a dihedral group, de- 
fined by the matrices: 



J(0) = 



V2 V2 

V2 \/2 
i 



(512) 
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(506) 



J4(0) 



-10 
0-10 
1 



(513) 



The particular matrices we have been considering above 
play the role of the phase gate, one of the possible equiv- 
alent C-NOT trit gates and the Hadamard gate in time 
optimal qutrit quantum computation. They are signed 
complex permutation matrices. The extremals of the uni- 
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tary generators are at the points {0, — , tt, — } which de- 
fine a base set of qutrit computational gates: 



0(0) 



^/2 ^ 



\/2 \/2 




Q{-. 



1 1 

? f 

V2 V2 





ie-' 
y/2 \/2 
\/2 V2 






-1 

(514) 



(515) 
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V2 V2 

i —i 



-i 
' 





1 

"71 71 



^2 





(520) 

This set of unitary transformations is necessary and suf- 
ficient to implement any quantum trinary computation. 
Some elementary operations include: 



^(0)^(|)-^(|)^(0) = 



£)(0)^(|)+^(|)^(0) 



i i 
i -i 




1 -1 
1 1 
2 



(521) 



(522) 



0(|) 



-ie^'e 

1 1 

V2 x/2 



V2 V2 







(516) 



It is possible to classify these matrices by their eigen- 
values and determinants into various point symmetry 
groups. 

In some sense the unitary matrices we are using play 
the role of "square roots" of unitary matrices themselves, 
which are in fact the dihedral group above. Being the 
complete set of positive and negative roots of the group 
expression, we receive twelve unique operators which may 
be readily used to generate a wide variety of useful oper- 
ations. 



J(0) 
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\/2 ^ 

^/2 ^/2 
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Jin) 
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- t 

V2 V2 

-i 



(517) 



Other operators will be explored in future works. By 
classification of the eigenvalues, these matrices may be 
defined as a certain types of point and continuous sym- 
metry groups in the complex plane, which are related to 
the roots of polynomial equations which we explored ear- 
lier. This method of qutrit computation and modelling 
of quantum states and control is efficient, in that it is 
time optimal; resource bounded, in that there is an ex- 
isting physical overhead on the total energy imparted to 
the system, reversible, in that it is unitary, and robust to 
error, as it is periodic. We could consider our SU(3) sys- 
tem to be in all ways equivalent to the physical sum of a 
qubit and an ancilla. In this sense, we are moving beyond 
qubit quantum computation into new realms of possibil- 
ity. With a qubit and an ancilla it is possible to achieve 
much more than with a single qubit, or any number of 
coupled qubits. Data correction requires ancillas, and 
the nature of quantum states is that all indistinguishable 
alternatives are to be summed over; this necessitates our 
examination of qutrits both as quantum computational 
objects in themselves, and as possible sources of error 
within a qubit calculation. This is the simplest example 
of a odd dimensional quantum state; in this particularly 
simple case we are fortunately able to derive rich symme- 
tries which are descriptive of many of the full properties 
of groups which do not emerge in lower dimensional cases. 
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XXV. SUBGROUPS OF SU(3) 

Explicitly writing out the diagonal generators of SU(3), 
we find: 



1 

7! 

("523) 

The diagonal subgroup has one component which com- 
mutes with every other member of the group, the other 
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only commutes with the rotation subgroup. We may 
therefore use the diagonal component to define the group 
structure, as it does not matter whether the commuting 
operator is taken as part of the constraint or Hamilto- 
nian. The matrix structure of the possible subdivisions 
is then: 



(524) 



vector relationship: 



d 
di 
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(531) 
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Our time optimal Hamiltonian is therefore a constant 
matrix, and may be exponentiated directly. Our wave- 
vector after some time is given by: 



|V(<)) = exp(-im) IV'(O)) 



(532) 



(525) 



Evaluating the quantum brachistochrone equation, and 
choosing the initial state |'!/'(0)) = [1,0,0,0]"^ we imme- 
diately find A;;; = — -^j/- Our final state is: 

|V'(r)) = i=(|00)-z|ll)) (533) 
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from which we obtain ~ 0. Our state evolves accord 


K* 













£2 


W3 


ing to: 



' i9 


/3 


o" 




7 





£1 






/3* 












7 


£2 


(527) 















. '^1 


£2 


-27 _ 







1 

71 



cos{2Xxt) 





(534) 



For any choice of Hamiltonian operator and associated 
constraint, we either find one or the other to be peri- 
odic (or constant). The equations of dynamic state are 
equivalent up to isometry with a permutation matrix. 



XXVI. GEODESIC ON SU(4) 

Let us examine a system where we can only apply 
global control pulses, without individual addressing, on 
two qubits. The Hamiltonian matrix is the Heisenberg 
model: 



H[t] = Xx^X ® (yx+ Xy&y ® &y + X^^Z ^ (7 z (528) 

In matrix form: 



and therefore our energy-time relationship is T 



XXVII. DIRAC EQUATION IN CO-ROTATING 
FRAME 



We take the more general ansatz for a Hamiltonian on 
SU(4): 
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-Xz X+ 
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; A± = A, ± Xy (529) 



H = 



Evaluating the determinant: 

det{H-Xi) = 
Obtaining eigenvalue equation: 

{-X' + a^+pl + \e\'f^Q 



The constraint may be expanded in the Hermitian basis: 

F[t] = {m{t) -5)®! + !® {n{t) -3) + ^ ^ijit)a., ® aj 

(530) 

Evaluating the quantum brachistochrone, we find the 



X = ±^a^ +pl + \e\' 
The isotropic constraint gives: 

Tr{ — ) = const. = a'^ + pi + \ef 



(535) 

(536) 

(537) 
(538) 

(539) 
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We therefore rescale the time parameter such that 
+ p1 + \£\^ = 1, so our original Hamiltonian matrix 
has the form of a 4-dimensional unit vector. The Hamil- 
tonian matrix is difficult to work with, in strict terms 
of evaluating the quantum brachistochrone equation; the 
solutions are not obvious, and the bilinearity of the differ- 
ential equations causes difficulty in finding an exact so- 
lution. We therefore must use a unitary transformation, 
equivalent to the Hadamard gate on the Pauli matrices, 
that acts on the 4-spinor. The necessary transformation 
is given by: 



V2 



1 1 

i -i 



1 



(7^) «) i 



WW'' = W'^W = 1 



W = W'' 



(540) 



(541) 



(542) 



Using the standard symmetry transformation, we obtain 
a new Hamiltonian matrix: 



The unitary operator which transforms the Hamiltonian 
is of the form Q(s) = P{s)W . This has the useful prop- 
erty: 



[F^P,Q(s)] =0 



(552) 



which enables the direct evaluation of the time evolution 
operator: 



C/(t,0) 



ie** 
ie-^* 



(553) 



Using this unitary to transform the initial condition of 
the Hamiltonian matrix, we find: 

ii{t) = U{t, 0)i7(0)?7l"(i, 0) = Ho + V[t, /3(0)] (554) 



H{t) 



.(0)1 



*it(o) 



e2"^(0) -a(0)l 



H{t + T) = H{t) 



(555) 



(556) 



Hy^, = WHW^ = 
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6t 6 



b = al + ifj ■ a 

^ ai — i/3 ■ a 
We then form matrices of eigenvectors: 
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The time evolution operator for the original system may 
be written as: 



U{t,0) = exp{-i / H{s)ds) 
Jo 

In the transformed reference frame this reads as: 



(550) 



det{H{t) - Ai) = det{H{Q) - Al) = 



(557) 



The matrix constraint for the time optimal control prob- 
lem can then be written as: 
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e2»*it 



6 

^2 



(558) 



where A 



a b 
c d 



and all the variables involved (other 



than the obvious time dependence) are constants. This 
Hamiltonian is the block-diagonal solution to the Klein- 
Gordon equation, in that we may write: 
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i 



(559) 



Obviously this is invariant under the time dependent 
transformation on the Hamiltonian. It is possible to map 
this solution unitarily to a number of other physical sys- 
tems which describe the motion of electrons and positrons 
in different basis sets. For completeness, these matrices 
are listed below. 
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U(t,0) = exp{~i Q\s)HwpQis)ds) (551) 
Jo 
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These transformation matrices mean that, given we have 
solved the specific problem of the Dirac matrix, that a 
large set of problems fit into the equivalence class, which 
saves a great deal of effort developing particular solu- 
tions. Note also that the Hamiltonian has a certain pe- 
riodicity, which can be written as Tmin x \\E\\ = tt, this 
relationship can be extended unitarily to all the other 
equivalent representations by isometry. 



XXVIII. DIMENSIONAL ARRAYS OF 
ARBITRARY SIZE 

Take an arbitrary tracefree matrix in SU(n), which we 
can express as the array: 



X = 



• ■ • a„_i 
an 

a* UJr. 



, ^ cjj = (569) 



This can be mapped invertibly and one-to-one to the 
vector 



a;„), (ai, ••• , a^), (aj, ••• , a^)] 

(570) 

which exists in the space [M"-i, C^ C''^* ] , and the 

dimension of the complex space is A; = X]""^ • O^'^ of 
the components of the real vector within the 
multivector is dependent; we retain the dummy variable 
to aid in the analysis. In any real calculation we 
immediately remove the redundant variable once the 
differential equations are calculated. If we choose the 
particular Hamiltonian matrix: 
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H = Y,ujn{t) \n) {n\ + (e,(t) |fc) (j| + h.c.) (572) 

Then the constraint is given by: 

F= J2 («.,,(*) N) 01+ h.c.) (573) 
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/32,„ 



/5n-l,n 
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We may write the quantum brachistochrone as the form: 
d 



Hence the Hamiltonian is constant: 







(575) 



(576) 



Note that, by induction, we have shown the existence 
of a geodesic on SU(n), having shown the existence of 
geodesies on SU(2), SU(2+1), SU(n) and SU(n+l). An- 
other family of infinite dimensional matrices is given by: 
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with associated constraint matrix F 
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(578) 

Again, using the quantum brachistochrone, we find that 
the Hamiltonian matrix, which defines a continuous sym- 
metry group, is constant. One more infinite dimensional 
system has system Hamiltonian and constraint: 
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F[t] 



Vl,2 



Vl.n 



77,i-i,„ 

Vn-l.n 
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j{t) = c<Ji+aj2 + .-.-+Wj+....+i^,i-i+w„ = (581) 



Obviously, by construction, the Hamiltonian commutes 
with the constraint, and is therefore constant. The diffi- 
culty with these problems on arbitrary dimensions is that 
there is no obvious way in which the Lagrange multipli- 
ers, representing total upper bounds on the energy, gener- 
alises to a matrix of infinite dimensions. One would hope 
that the sums would go over to integrals in the correct 
fashion, but this remains to shown. Of course it is possi- 
ble to create these type of general Hamiltonian-constraint 
systems with the geodesic property, but it would be de- 
sirable to have expressions for field densities for matrices 
of infinite dimension, which represent the continuous de- 
grees of freedom. 

In conclusion, we have shown within this paper how one 
might go about calculating these more diflacult 
examples on finite dimensional systems, as well as 
demonstrating the links between the matrix and vector 
dynamic calculus. This paper has demonstrated that 
the use of finite matrix methods and algebraic 
techniques may yield dividends, even when the situation 
is graphically complex. We have used the quantum 
brachistochrone equation to derive a number of new 
results on SU(3); these methods may be extended to 
higher-dimensional matrix groups in order to find 
time-optimal flows of quantum states. The results 
derived are intriguing, especially the links between the 
various branches of mathematics. Our principal results 
may be applied within the laboratory to achieve fast 
data transfer; as the method of calculation is not reliant 
on adiabatic transfer or the rotating wave 
approximation it is theoretically the most general and 
efficient scheme that may be applied. 
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